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ABSTRACT 

The 30 Myr old A3-type star HD 21997 is one of the two known debris dust disks having a measurable amount 
of cold molecular gas. With the goal of understanding the physical state, origin, and evolution of the gas in young 
debris disks, we obtained CO line observations with the Atacama Large Millimeter/submillimeter Array (ALMA). 
Here, we report on the detection of 12 CO and 13 CO in the J = 2-1 and J = 3-2 transitions and C 18 0 in the J — 
2-1 line. The gas exhibits a Keplerian velocity curve, one of the few direct measurements of Keplerian rotation in 
young debris disks. The measured CO brightness distribution could be reproduced by a simple star+disk system, 
whose parameters are r m < 26 AU, r out = 138 ± 20 AU, M* = 1.87 q 2 ^©. and ' = 3216 ± 3: 1 . The total CO 
mass, as calculated from the optically thin C 18 0 line, is about (4-8) x 1 0 2 M®, while the CO line ratios suggest a 
radiation temperature on the order of 6-9 K. Comparing our results with those obtained for the dust component of 
the HD 21997 disk from ALMA continuum observations by Moor et ah, we conclude that comparable amounts of 
CO gas and dust are present in the disk. Interestingly, the gas and dust in the HD 21997 system are not colocated, 
indicating a dust-free inner gas disk within 55 AU of the star. We explore two possible scenarios for the origin of 
the gas. A secondary origin, which involves gas production from colliding or active planetesimals, would require 
unreasonably high gas production rates and would not explain why the gas and dust are not colocated. We propose 
that HD 21997 is a hybrid system where secondary debris dust and primordial gas coexist. HD 21997, whose age 
exceeds both the model predictions for disk clearing and the ages of the oldest T Tauri-like or transitional gas disks 
in the literature, may be a key object linking the primordial and the debris phases of disk evolution. 
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1. INTRODUCTION 

Nearly all young stars harbor circumstellar disks whose ther- 
mal emission produces a strong infrared excess with fractional 
luminosities of L^/L ^ 0.1. At the early stage of their evolution, 
the masses of these disks (typically a few percent of the mass of 
the central star) are dominated by gas with only a few percent 
of the mass in small, submicron-sized dust grains. Both the dust 
and the gas are of primordial origin. As the disk evolves, dust 
settles in the midplane and the grains eventually form larger bod- 
ies, planetesimals, and planets. The gas is removed by viscous 
accretion (e.g., Lynden-Bell & Pringle 1974), by photoevapo- 
ration (e.g., Alexander 2008), or by planet formation. Current 
observational results imply that the primordial gas is mostly de- 
pleted at ages of <10 Myr (Pascucci et al. 2006; Fedele et al. 
2010). 

A detectable amount of infrared excess is also present in many 
older, main-sequence stars, with typical fractional luminosities 
of La/L ^ 10~ 3 (e.g., Roccatagliata et al. 2009). These debris 
disks are fundamentally different from primordial disks. Their 
masses, as inferred from the emission of small dust grains, 
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are usually below 1 M® and they are practically gas-free. 
Without the stabilizing effect of surrounding gas, the lifetime 
of individual dust grains in debris disks is very short due to 
removal by dynamical interactions with stellar radiation. Thus, 
dust needs to be continuously replenished by collisions and/or 
evaporation of previously formed planetesimals (e.g., Wyatt 
2008). The same processes would in principle produce gas as 
well, via sublimation of planetesimals (Lagrange et al. 1998), 
photodesorption from dust grains (Grigorieva et al. 2007), 
vaporization of colliding dust particles (Czechowski & Mann 
2007), or collision of comets or icy planetesimals (Zuckerman 
& Song 2012). Thus, the secondary gas that may be present in a 
debris disk is dominated by CO and HiO (Mumma & Charnley 
2011) and only a small amount of Hj, mainly originating from 
the dissociation of HiO, is expected. 

In theory, it is possible that systems in transition between 
the primordial and the debris state possess hybrid disks, i.e., 
primordial gas is accompanied by secondary dust. It is also 
possible that the outer disk is still primordial, while the inner 
disk is already composed of secondary material (Wyatt 2008; 
Krivov et al. 2009). So far, only six debris disks with detectable 
gas components are known. The edge-on orientation of the disks 
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Table 1 

Rest Frequencies, Observing Setups (See Text for Details), Beam Sizes and Position Angles, Noise Levels, and Line Fluxes for the CO Observations of FID 21997 


Line 

Frequency 

(GHz) 

Setup 

Beam Size 
("x") 

Beam RA. 

(°) 

Noise 

(mJy beam -1 channel -1 ) 

Peak Flux 

(Jy) 

Total Flux 
(Jy km s -1 ) 

Total Flux 
(W m~ 2 ) 

12 CO(3-2) 

345.796 

Bl 

1.19 x 1.47 

117.5 

7.9 

0.615 

2.52 ± 0.27 

(2.91 ± 0.31) x 10~ 20 

12 CO(2-l) 

230.538 

A2 

1.14 x 1.21 

16.4 

8.6 

0.566 

2.17 ±0.23 

(1.67 ± 0.18) x 10~ 20 

13 CO(3-2) 

330.588 

B1.B2 

1.20 x 1.38 

-49.9 

5.7 

0.387 

1.35 ± 0.15 

(1.49 ±0.17) x 10~ 20 

13 CO(2-l ) 

220.399 

Al 

1.19 x 1.92 

19.0 

7.6 

0.238 

0.82 ±0.10 

(6.04 ± 0.74) x 10“ 21 

C 18 0(2-1) 

219.560 

Al 

1.21 x 1.92 

18.6 

9.0 

0.129 

0.36 ± 0.06 

(2.63 ± 0.44) x 10“ 21 


around f J > Pic and HD 32297 allow the detection of a very 
small amount of circumstellar gas based on the presence of 
absorption lines (Roberge et al. 2000; Redfield 2007). Another 
member of the /3 Pic moving group, HD 172555, also contains 
some gas, as evidenced by [O i] emission at 63 /im (Riviere - 
Marichalar et al. 2012). Recently, the [C n] line at 158 )i m was 
detected in ~30 Myr old HD 32297 (Donaldson et al. 2013). 
A substantial amount of molecular gas has been detected at 
millimeter wavelengths in the debris disk around the young 
main-sequence stars 49 Ceti (Hughes et al. 2008) and HD 21997 
(Moor et al. 2011). The ages of these systems are between 10 
and 40 Myr, which partly overlaps with the transition period 
from the primordial to the debris phases. Although HD 21997 
and 49 Cet are the oldest members of the gaseous debris disk 
sample, they also contain the largest amount of cold gas. Thus, 
it is worth considering whether in their case the detected gas is 
of primordial or secondary origin. 

In this paper, we focus on HD 21997, which is an 
A3-type star at a distance of 72 pc (based on Hipparcos par- 
allax; van Leeuwen 2007), belonging to the well-dated 30 Myr 
old Columba moving group (Moor et al. 2006; Torres et al. 
2008). In an earlier study, we detected molecular gas in the 
J = 3-2 and J = 2-1 transitions (CO(3-2) and CO(2-l)) with 
the APEX telescope (Moor et al. 2011). Motivated to under- 
stand the origin and evolutionary status of the gas component 
in this system, we obtained (sub)millimeter interferometric 
continuum and CO line observations with the Atacama Large 
Millimeter/ submillimeter Array (ALMA). Our aim was to spa- 
tially resolve the disk in order to study the relative location of 
the gas and dust components, determine the distribution of the 
different CO isotopologues, precisely measure the gas mass, 
and analyze the gas kinematics. In this paper, we present our 
results on the gas content of the disk, while the dust contin- 
uum observations are analyzed in Moor et al. 2013 (hereafter, 
Paper I). 

2. OBSERVATIONS AND DATA REDUCTION 

We observed HD 21997 with ALMA in Cycle 0 using a 
compact configuration (PI: A. Kospal). The quasar J0403— 360 
served as a bandpass and gain calibrator, while Callisto was 
used to set the absolute amplitude scale using the CASA 
Butler-JPL-Horizons 2010 model. We obtained both line and 
continuum data. Here, we focus on the line observations, while 
the continuum measurements are presented in Paper I. 

In Band 6, we targeted the (2-1) transitions of 12 CO, 13 CO, 
and C ls O. We used the frequency division correlator mode 
(FDM). Each of the four simultaneously observed spectral 
windows offered 3840 channels with a channel separation of 
122 kHz, resulting in a bandwidth of 469 MHz per window. 
The spectral resolution corresponds to a velocity resolution 
of 0.33 km s -1 at 230 GHz. Due to restrictions in the setup 
of the local oscillators, the three CO lines could not be 


observed simultaneously; therefore, two different correlator 
configurations were prepared. One setup (Al) was configured to 
observe °CO(2-l) and C ls O(2-l) and it was executed on 201 1 
November 29 using 14 antennas, with baselines ranging from 9.7 
to 148 kA. The other setup (A2) was tuned to have 12 CO(2-l) 
in one of the spectral windows and it was observed on 2011 
December 31, using 16 antennas, with baselines between 9.6 
and 196 kA. Each of these two setups had an on-source time of 
49.4 minutes. 

In Band 7, we targeted the (3-2) transitions of 12 CO and 
13 CO, as well as the (7-6) transition of CS. The correlator 
was set up in FDM with 3840 channels, a channel separation 
of 244 kHz, and a total bandwidth of 938 MHz. The spectral 
resolution corresponds to a velocity resolution of 0.44 km s“ 1 
at 345 GHz. Again, two correlator configurations were used. 
In one setup (Bl), 12 CO(3-2) and 13 CO(3-2) were targeted. 
Observations were done on 2011 November 3 and 4 using 14 
and 15 antennas, respectively, with baselines ranging from 14 to 
154 kA. The other setup (B2) covered the CS(7-6) line and was 
observed on 2011 November 3 with 14 antennas and baselines 
between 14 and 153 kA. The rest frequency of CS(7-6) also 
allowed us to obtain additional coverage of 13 CO(3-2) in this 
configuration. Both configurations had an on-source time of 
49.4 minutes, resulting in ~ 1 00 minutes of observing time for 
13 CO(3-2). 

We subtracted the continuum emission in the spectral cubes 
using the CASA task uvcontsub and individually cleaned the 
spectral region around each molecular line using the CASA 
task clean with robust weighting. A summary of the observing 
setup, beam size and position angle (PA.) for each line, and 
the rms noise measured in the channel maps is presented in 
Table 1. 

3. RESULTS 

3.1. Integrated Line Profiles 

We detected all targeted CO lines, i.e., 12 CO(2-l), 12 CO(3-2), 
13 CO(2-l), 13 CO(3-2), and C ls O(2-l). Neither the CS(7-6) 
line nor any other lines are detected in our ALMA data cubes. 
In order to obtain integrated profiles for the CO lines, we added 
all the flux within a radius of 2" 8 of the stellar position (the 
smallest aperture that contains the whole flux of the object) for 
each velocity channel. The resulting line profiles are plotted in 
Figure 1 (left panel). In line with our earlier APEX observations 
(Moor et al. 2011), the lines are all double-peaked, typical for 
emission arising from rotating material. The intensities of the 
two 12 CO lines are also consistent within the uncertainties with 
our earlier single-dish APEX results, indicating that there is no 
significant flux loss due to the interferometer filtering out large- 
scale structures. The peak and central velocities are the same 
for all lines. Within the uncertainties, the profiles are symmetric 
around a systemic velocity of u sys = 1.29 km s" 1 (measured 
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Figure 1. CO spectra of HD 21997 observed with ALMA. The left panel shows integrated flux densities in Jy, while the right panel displays spectra normalized to 
their respective line areas, to better show the similarities in the line profiles. In the right panel, an additional thick gray line shows the best-fit model described in 
Section 3.4. 

(A color version of this figure is available in the online journal.) 
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Figure 2. Left panel: zeroth moment (with contours) and first moment (with color scale) of the 2 CO (2-1) emission. The major and minor axes are plotted with 
straight black lines. The gray ellipse in the lower left corner indicates the beam. Right panel: position-velocity diagram of 12 CO (2-1) emission along the major axis. 
Overplotted in white is a Keplerian rotation curve for M, = 1.8 Mq and i = 32:6 (see Section 3.4). 

(A color version of this figure is available in the online journal.) 


in the local standard of rest (LSR) system), which agrees with 
the radial velocity of the central star. The intensity-weighted 
average velocity for all lines is the same (1.29 ± 0.05 km s -1 ) 
as well. The symmetric profiles suggest that the rotating material 
is distributed axisymmetrically around the central star. 

Table 1 presents the line fluxes integrated from —4 km s _ 1 to 
6.6 km s -1 , a velocity interval that covers the whole line, as well 
as the peak fluxes. The strongest CO line is the 3-2 transition 
of 12 CO, followed by the 2-1 line of the same isotopologue, 
then the 3-2 and 2-1 transitions of 13 CO, while the 2-1 line 
of C 18 0 is the faintest among our detections. If we normalize 
the line profiles (Figure 1, right panel), we find that the lines 
have very similar profiles. We checked which velocity channels 
contain significant emission above the 3a level and found that 
these span a range of ±4.3 km s -1 around the systemic velocity. 

3.2. Spatially Resolved CO Emission 

The disk was spatially resolved in all CO transitions. We 
calculated zeroth and first moment maps for each data cube 
(for an example, see Figure 2, left panel). Figure 2 shows that 
the intensity profile peaks at the stellar position and smoothly 
decreases with radial distance. The velocity maps indicate that 
to the north of the star, material is approaching, while in the 
south, material is receding, with respect to u sys . The same 
characteristics can be seen for the other lines as well. Assuming 
rotating disk kinematics, we can determine the rotation axis by 
fitting a line to those pixels in the first moment map where the 


velocity equals v sys . This axis also marks the minor axis of the 
inclined disk image. We determined the P. A. of the major axis on 
each image and obtained an average of 22? 6 ± 0?5. The extent 
of the emission along the minor and major axes defined by the 
3a contour in Figure 2 (left panel) is 3'.'5 x 4'.'0. 

Figure 2 (right panel) shows a position velocity diagram 
measured along the major axis of the disk marked in Figure 2 
(left panel). Most of the emission comes from two well- 
separated regions. The highest velocities are observed closest 
to the stellar position and velocity gradually decreases with 
increasing distance, a pattern characteristic of Keplerian-like 
disk rotation. In Figure 3 (left panel), we plot the observed CO 
emission in nine different velocity channels between —3.2 and 
3.2 km s -1 relative to i> sys , in steps of 0.8 km s _1 . The plots show 
that the emission in the most extreme velocity channels is rather 
compact and is located close to the stellar position. Toward 
lower velocities, the peak of the emission appears farther from 
the stellar position. For velocities close to the systemic velocity, 
the emission is double-peaked. The morphology of the images is 
very similar on either side of the systemic velocity, apart from a 
mirroring. It is evident from the plots that the different CO lines 
exhibit the same emission pattern, allowing for the differences 
in the beam shape/size and the signal-to-noise ratio. 

3.3. Temperature and Mass of the CO Gas 

Despite their very different abundances, the fact that the 
intensities of the three different CO isotopologues are of the 
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Figure 3. Left: channel maps of CO observations of HD 21997. The velocities indicated in the upper left comers are relative to the systemic velocity. The solid 
contours mark the 3 a, 6a, 9a, 12 a, etc... levels, while the dashed contours indicate the — 3cr level. Right: observed and modeled CO channel maps for the 12 CO (2-1) 
line and residuals. Details of the model are described in Section 3.4. The contours are the same as in the left panel. 

(A color version of this figure is available in the online journal.) 


same order of magnitude already indicates that 12 CO is probably 
optically thick. We observed the ( 2 - 1 ) transition for three 
different CO isotopologues. If we denote the optical depths 
of 12 CO, 13 CO, and C 18 0 by i\2, T13, and ns, respectively, then 
the ratio of the 12 CO to the C 18 0 line can be approximated 
as (1 — e _Tl2 )/(l — e ~ z,s ) and a similar formula holds for the 
13 CO to the 18 CO line ratio. Assuming that the optical depths of 
the different isotopologues follow the same proportions as the 


abundance ratios typical of local interstellar matter (Wilson & 
Rood 1994), then t\2 = 560tig and T13 = 7.4ti8. Using these 
numbers, we obtained ris = 0.2 from the 12 CO to the C 18 0 line 
ratio and x\% = 0.6 from the 13 CO to the C ls O line ratio. This 
result means that T13 is in the 1.5-4. 5 range, while ti 2 is between 
about 100 and 300. Using comet-like isotopic ratios would yield 
similar results (Eberhardt et al. 1995; Bockelee-Morvan et al. 
2012; Rousselot et al. 2012). According to Visser et al. (2009), 
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and references therein, the 12 C/ 13 C ratio may be a factor of two 
smaller or larger, while the 16 0/ 18 0 ratio may be three to five 
times larger than the local elemental ratios. This fact means that 
we may overestimate r® and ri 2 by less than a factor of two, or 
we may underestimate them by as much as a factor of five. 

Since the 12 CO lines are optically thick, the ALMA 
12 CO(2-l) and 12 CO(3-2) maps can be used to calculate the 
radiation temperature of the gas by comparing the observed 
surface brightness with a blackbody. For optically thick lines, 
the radiation temperature equals the excitation temperature. 
Because the intrinsic line width and the significant macroscopic 
movements of the gas as it orbits the central star decrease the 
optical depth at a certain frequency, we always selected the 
frequency channel with the maximal intensity for each map 
pixel. This way, for each map pixel, we can estimate the tem- 
perature at the frequency (or velocity) where the optical depth 
is the highest. For 12 CO, we obtained excitation temperatures 
up to 9 K. The temperature can also be estimated from the 
ratio of the optically thick 12 CO(3-2) and 12 CO(2-l ) lines. In the 
Rayleigh-Jeans approximation, the ratio is expected to be the ra- 
tio of the squares of the line frequencies, i.e., about 2.25. Instead, 
using the total flux values from Table 1, we obtain 1.16. The 
low value suggests that the temperature is very low and that 
the Rayleigh-Jeans approximation is not valid. Indeed, using the 
Planck function, it turns out that these line ratios correspond to 
5.6 K. 

In Moor et al. (201 1), given a lack of more information, we 
assumed that 12 CO was optically thin and estimated a total CO 
mass of 3.5 x 1 0 4 M®. Our ALMA observations made it evident 
that l2 CO was optically thick and that C 18 0 should be used for 
mass estimation. The calculations depend on the temperature 
of the gas, for which we used different values between 5.6 K 
(the excitation temperature obtained from the 12 CO (3-2) to 
(2-1) line ratio) and 100 K (the hottest dust temperatures at 
the inner edge of the disk in the model of Paper I). Higher 
temperatures mean higher emission, but also mean that these 
low -J transitions are less excited. The net result of these two 
effects is that within the assumed temperature range, the total 
CO mass needed to produce the observed emission changes by 
only a factor of two. Thus, the CO mass is well constrained, even 
without knowing the precise temperature. The resulting value, 
(4-8) x 1 0“ 2 ;V/q, is about two orders of magnitude higher than 
previously believed. Considering the uncertainties in the 12 CO, 
13 CO, and C 18 0 isotopic ratios discussed above, the total CO 
mass may actually be a factor of two smaller or a factor of five 
higher than this value. Nevertheless, the precise value does not 
change the fact that HD 21997 has an unusually gas-rich disk. 

If the gas has a second-generation origin, the most common 
species, based on cometary composition, are H 2 O, CO, and 
CO 2 in comparable abundances (Mumma & Charnley 2011 and 
references therein). Thus, the CO mass gives a good estimate for 
the order of magnitude of the total gas mass in the disk. However, 
if the gas is primordial, then not only CO but also H 2 gas is 
present in the disk and the total gas mass is significantly higher. 
Taking a canonical CO/H 2 abundance ratio of 10 4 , the total gas 
mass is on the order of 26-60 M®. If such a large amount of H 2 
is indeed present, and it is warm enough, the gas should display 
rotational lines at mid-infrared wavelengths. We checked the 
Spitzer /Infrared Spectrograph spectrum of HD 21997 (Paper I) 
for the 5(0), 5(1), and 5(2) lines, but found that they are not 
present. Using a 3 a upper limit for the flux of the line expected 
to be the brightest, 5(0), we calculated an upper limit of 35 M® 
for the total gas mass for 100 K gas and 3400 M® for 50 K gas. 


This result means that even if H 2 is present, it should be quite 
cold, just like the CO gas. 

3.4. Distribution of the CO Gas 

To determine the basic parameters of the HD 21997 gas disk, 
we fit the spatial and velocity distributions of the observed CO 
emission with a simple disk geometry combined with a Keple- 
rian velocity profile. Our simple model has five parameters: the 
stellar mass (M*), the inner and outer disk radii (r m , r out ), the 
inclination of the disk (i, where i — 0° corresponds to face on), 
and the exponent (p ) of the radial brightness profile, assumed 
to be a power law (/ (r) oc r p ). For the velocity held of the gas, 
we adopted Keplerian rotation around the central star, expected 
to be true for low-mass, non-self-gravitating disks (e.g., Dutrey 
et al. 2007). We used an intrinsic line width of 0.1 km s , corre- 
sponding to the Doppler broadening of a line arising from 10 K 
gas. We calculated the models for the same velocity channels as 
the actual observations and applied two-channel wide Hanning 
smoothing to account for the instrumental velocity resolution. 
In the spatial direction, we first calculated the models with a six 
times oversampling, convolved the models with the beam of the 
actual observations, and rebinned them to the same pixel scale 
as the observations. We decided to model the 12 CO (2-1) line, 
which is one of the brightest we observed and has the smallest 
and most circular beam. 

To explore a sufficiently large parameter range, we changed 
M* between 0.9 M® and 3.0 M® in steps of 0.3 M®, r m be- 
tween 5 AU and 60 AU in steps of 5 AU, r out between 95 AU 
and 175 AU in steps of 5 AU, i between 25° and 39° in steps 
of 2°, and p between 0 and —1.75 in steps of 0.25. In total, 
we calculated 104,448 models. In order to compare the mod- 
els with the observations, we normalized the models so that 
the total intensity equaled the observed total intensity. Then, 
we computed the / 2 for each channel map. The final / 2 of a 
certain model was debited as the maximum of the x 2 of all 
channels. Finally, we calculated the Bayesian probability for 
each model as exp(— x 2 /2) and normalized them so that the 
sum of the probabilities of all models is one. Figure 4 shows the 
one-dimensional (ID) probability distributions as a function of 
each parameter and the two-dimensional (2D) probability dis- 
tributions as a function of two different parameters. These are 
marginal distributions, i.e., for example, in the distribution as a 
function of r m , all the other parameters (r out , M„, i, and p ) were 
marginalized out. Similarly, for example, in the case of the 2D 
distribution for r m and r out , M*, i, and p were marginalized out 
(for more information on Bayesian inference and marginaliza- 
tion, see, e.g., Loredo 1992). The value where the ID probabil- 
ity distribution is maximal formally gives the best estimate for 
that parameter. Formal la uncertainties can also be estimated 
from these graphs by integrating the probability distributions 
on each side of the maximum until one obtains 0.68. However, 
it is evident from Figure 4 that not all parameters are equally 
well determined. The probability distribution as a function of 
n n is quite flat for short inner radii, indicating that we can only 
give an upper limit on r ln . The outer radius r out and stellar 
mass M* are well determined, while the power-law exponent 
of the brightness profile, p, shows a very slowly changing and 
wide probability distribution, which practically makes this pa- 
rameter ill-determined. The inclination and the stellar mass are 
degenerate, thus we did not determine the inclination from the 
Bayesian probabilities, but calculated it directly from the best-fit 
stellar mass. 
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Figure 4. Marginal probability distribution for the different model parameters. Details of the model are described in Section 3.4. 


Table 2 

Parameters of the HD 21997 System and Its Gas Disk 


Distance 

d (pc) 

71.95 

Systemic velocity 

«LSR (km s _I ) 

1.29 

Disk orientation 

pa. n 

22.6 ±0.5 

Disk inner radius 

On (AU) 

<26 

Disk outer radius 

0,ut (AU) 

138 ± 20 

Stellar mass 

M* (Mq) 

i o+0.5 
A ’°— 0.2 

Disk inclination 

<(°) 

32.6 ±3.1 

Disk brightness exponent 

P 

-1.1 ± 1.4 


Table 2 shows our best estimates for the disk parameters 
and their la uncertainties, along with a la upper limit for r m . 
Figure 3 (right panel) shows side-by-side the channel maps for 
the 12 CO (2-1) observations, the model with the best parameters 
(and r m — 20 AU), as well as the residuals. Moreover, we plotted 
the integrated line profile from the model in Figure 1. These 
figures demonstrate that the best model we found indeed fits 
the channel maps and line profile very well. We convolved this 
model with a beam appropriate for our other CO observations 
and scaled the models so that the integrated line intensity equals 
the observed intensities. We then plotted the channel maps and 
residuals similarly to those in Figure 3 (right panel) and found 
that the observations and the model match each other well. 
This fact demonstrates that a single-disk model with one set of 
parameters reproduces all five CO lines. 

In the following, we check whether or not the fitted parameters 
are consistent with our a priori expectations. The major axis of 


the 12 CO(2-l) image is 4"0 (Section 3.2), which, deconvolved 
with the beam size of 1"21 (Table 1), gives an outer radius of 
137 AU, in agreement with our fitted value. We can estimate 
the inclination from the ratio of the minor and major axes 
of the 12 CO(2-l) image, which gives 29° (cf. 32?6 ± 3?1 in 
Table 2). A similar calculation was done from the minor and 
major axes of the ellipse fit to the ALMA dust continuum image 
(Paper I), giving 32°9 ± 2°6. Another argument in favor of a 
~3(T inclination is that from the measured v sin i = 70 km s' 1 
(Royer et al. 2007), we obtain an equatorial rotational velocity 
of 140 km s -1 , typical for A-type stars with a mass of 
about 1.8 M q (Zorec & Royer 2012). Finally, the dynamical 
stellar mass of 1.8 M e agrees well with the mass of 1.85 M Q 
obtained from evolutionary tracks by Moor et al. (201 1). These 
comparisons demonstrate that our simple fitting procedure 
provided reasonable results for the basic properties of the CO 
gas distribution. A more detailed physical model, taking into 
account the vertical disk structure and full radiative transfer, is 
postponed to our forthcoming paper. 

4. DISCUSSION 

Our modeling in Section 3.4 suggests that the observed CO 
emission can be well described by a gas disk in Keplerian 
rotation around a 1.8 Mq central star. For the inner radius 
of the gas disk, we obtained a la upper limit of 26 AU. In 
Paper I, we found that the dust disk starts at a radius of »55 AU. 
Consequently, we can claim with a 99% confidence level that 
the inner radius of the gas is closer to the central star than that 
of the dust disk. This result is in accordance with the strikingly 
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different appearance of the ring-like dust continuum image and 
the centrally peaked CO zeroth moment maps. There is a part 
in the disk where the gas and the dust are not co-located and 
this “inner disk” is practically dust-free. The total CO gas mass 
is about (4-8) x 1 O' 2 M®. In Paper I, we modeled the dust 
emission of HD 21997 and found that the dust mass within about 
150 AU is 0.09 M®, resulting in a CO-to-dust ratio of 0.4-0. 9, 
i.e., roughly the same mass of dust and CO gas is present in the 
disk. This value may actually be slightly lower if we take into 
account that a fraction of the gas is located in the inner dust- 
free area. While the dust properties of the HD 21997 system 
are characteristic of a typical debris disk (Paper I), there is an 
unexpectedly large amount of gas in the disk (for comparison, 
the CO-to-dust ratio in /) Pic is less than 2 x 10 -4 ; Dent et al. 
1995; Nilsson et al. 2009). In the following, we discuss different 
possibilities for the origin, physical properties, and possible 
evolution of the gas component. 

4.1. Debris Disk with Pure Secondary Gas 

In a debris disk, where neither dust grains nor H 2 provide 
enough shielding for CO, the main factor determining the CO 
lifetime is its self-shielding against the stellar UV radiation field 
and the interstellar radiation field. To estimate the CO lifetime 
in the HD 21997 disk, we distributed the measured CO mass 
between 20 AU and 138 AU, adopted a power-law radial density 
distribution (n oc r - “, with a — 1.5, 2.0, and 2.5), and assumed 
three different vertical scale heights (H/r = 0.05, 0.1, and 
0.2). We computed the radial and vertical column density of 
CO in the disk and estimated the shielding factors for 12 CO 
and C 18 0 from the tabulated values of Visser et al. (2009). We 
neglected the shielding from the dust grains due to their low 
column density. We took into account both stellar UV flux and 
the interstellar radiation field, the same way as in Moor et al. 
(2011). The resulting lifetimes are below 30,000 yr for 12 CO 
and below 6000 yr for C 18 0. The Keplerian rotation of the gas 
and resulting Doppler shift of the line would further decrease 
the optical depth and reduce the lifetime. Thus, a replenishment 
of CO molecules is needed. 

Assuming an equilibrium between the dissociation and pro- 
duction of CO molecules and dividing the measured gas mass by 
the maximum lifetime computed for 12 CO, we obtained a lower 
limit for the gas production rate of about 10 19 kg yr " 1 (approx. 
1 0 9 M® yr -1 ). We emphasize that this production rate is a very 
optimistic lower limit, because for a significant part of the disk, 
CO lifetimes are lower than the value quoted above. Typical 
gas production rates in solar system comets are on the order of 
10 9 -10 10 kg yr -1 (Wyckoff 1982). Most of this gas, however, is 
H 2 O, and CO production is lower by a factor of 5-100 (Jewitt 
et al. 2007). This fact means that to keep up the observed CO 
mass in HD 21997, the continuous gas production of 10 10 — 10 1 1 
comets is needed. Alternatively, considering a larger comet like 
Hale-Bopp (total mass of 1.3 x 10 16 kg; Weissman 2007) and 
assuming a ~10% CO content (Sykes et al. 1986; Mumma & 
Charnley 2011), at least 6000 Hale-Bopp-like comets need to be 
completely destroyed every year. Unless one assumes a recent 
transient event in which a large amount of CO was produced, 
the calculated very high replenishment rate argues against the 
secondary origin of the CO gas in the HD 21997 system. 

The destruction of planetesimals may result in both debris 
dust and secondary gas (Section 1). If the gas and dust indeed 
have a common origin in the HD 21997 system, one would 
expect that they are colocated, unless some physical process 
separates gas and dust and moves the gas inward or the dust 


outward from where they were produced. One possibility is that 
the planetesimal ring starts about 55 AU from the star, as the 
location of the dust disk implies. Viscous accretion can transport 
material into the inner part, but it would act both on gas and the 
small dust particles. The other possibility would be to assume 
that the gas production via the destruction of comets occurs in 
the warm inner disk. However, comets would also produce dust, 
which was not observed. Thus, the existence of the dust-free gas 
disk cannot be explained in the secondary scenario. 

4.2. Primordial Gas in a Hybrid Disk 

The other possible scenario for the origin of the gas in the 
HD 21997 system is that the CO molecules are remnants of 
the primordial protoplanetary disk. Their survival would require 
very efficient shielding over the 30 Myr long history of the 
system. Attributing this shielding to the presence of H 2 gas 
(with a CO/H 2 abundance ratio of 10 -4 ) and distributing the H 2 
in the same simple disk geometry as we assumed for the CO 
in Section 4.1, we found CO lifetimes typically two orders of 
magnitude longer than without the presence of H 2 , long enough 
that some primordial gas can still be present in the system. 
There is, however, one important drawback of the primordial 
gas scenario. The calculation above shows that the H 2 number 
density in the midplane is high enough 0i H , > 10 6 cm 3 ) for the 
CO gas to be collisionally excited and be in local thermodynamic 
equilibrium (LTE), using the critical densities from Kamp & van 
Zadelhoff (2001). In this case, the excitation temperature equals 
the kinetic temperature. Our CO observations indicate very low 
excitation temperatures, much lower than the dust temperature 
(64 K, based on spectral energy distribution fitting done by 
Moor et al. 201 1). It is an open question what could cause such 
a difference. The relatively warm dust temperature explains why 
the CO does not freeze out onto the dust grains. 

The oldest phases in primordial disk evolution are represented 
by transitional disks. From a morphological point of view, 
HD 21997 might give the impression of a transitional disk 
where gas exists within the dust-free hole (for such examples, 
see Brown et al. 2009 or Casassus et al. 2013). However, 
we argue that HD 21997 is not a transitional disk because 
our results indicate that its dust content is debris-like. This 
conclusion is supported by the following. (1) In transitional 
disks, the gas-to-dust ratio is typically below the interstellar 
value of 100 (Keane et al. 2013), while in HD 21997, using the 
total gas mass including H 2 , this ratio would be on the order 
of 300-700. This result suggests that HD 21997 is unusually 
dust-poor, possibly because the majority of the dust is already 
locked up in planetesimals. (2) The dust mass is only 0.1 M®, 
significantly lower than typical for transitional disks. (3) There 
are no spectral signatures of small particles, either silicates or 
polycyclic aromatic hydrocarbons. (4) The limited lifetime of 
dust particles even in the presence of gas requires replenishment 
(Paper I). 

Based on these arguments, it is possible that we detected 
a hybrid system where the gas is primordial while the dust 
is already secondary. This result would also explain why the 
dust and gas are not co-located. Because the dust component 
of HD 21997 consists of large grains (>6/xm; Paper I), the 
grains trace the location of the planetesimal belt starting at about 
55 AU. The gas, however, is the remnant of the primordial 
circumstellar material and thus may fill the full disk. It is 
interesting to speculate about the further evolution of the gas 
component in this system. It may be a stable structure, or 
it may be in the phase of final disappearance. We note that 
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there are still unresolved issues with the hybrid disk scenario. 
What happened to the primordial dust in the central dust-free 
inner disk? Planetesimals already formed in the outer disk, 
as evidenced by the debris dust, but why are they missing in 
the inner disk where dynamical timescales are shorter? Is it 
possible that solids rapidly accreted to planet masses without 
affecting the gas component (by forming super-Earths rather 
than Jupiters)? The future discovery of similar hybrid systems 
may help to answer these questions. Some of the answers may 
also help to better understand the evolution of some transitional 
disks, where the central clearing may contain tenuous secondary 
dust, while the outer disks consist of primordial gas and 
dust. 

5. SUMMARY AND CONCLUSIONS 

We presented ALMA CO line observations of the disk around 
the 30 Myr old star HD 21997. We detected the (2-1) and (3-2) 
lines of 12 CO and 13 CO and the (2-1) line of C 18 0. The line 
profiles and channel maps show that the gas is in Keplerian 
rotation around the central star. We calculated a grid of simple 
gas disk models with varying inner and outer radii, inclinations, 
stellar masses, and radial brightness profiles of CO emission. 
Using Bayesian probability analysis, we found that the best 
model gives r m < 26 AU, r ou t = 138 AU, M* = 1.8 M®, and 
i = 32:6, while the brightness profile is undetermined. The CO 
line ratios and intensities suggest that the gas temperature is 
very low (on the order of 6-9 K). The total CO mass in the 
disk, as calculated from the optically thin C 18 0 line, is about 
(4-8) x 10” 2 M®. Comparing our results with those obtained for 
the dust component from the ALMA continuum observations in 
Paper I, we concluded that, in terms of mass, a similar amount 
of CO gas and dust is present in the disk. Interestingly, the gas 
and dust in the HD 21997 system are not colocated: there is 
an inner, dust-free gas disk extending from closer than about 
26 AU to about 55 AU from the star. 

We discussed two possible scenarios for the origin of the 
gas in the HD 21997 disk. First, we explored the possibility 
that the gas is of secondary origin, mainly CO produced by 
planetesimals. In this case, the sub-critical gas densities would 
lead to non-LTE conditions, which might explain the low 
excitation temperatures. However, the short CO lifetimes and the 
necessary high CO production rates exclude this scenario. The 
other possibility is that the gas is primordial. In this case, there 
is a large amount of Hi gas is the disk, leading to significantly 
longer CO lifetimes. A primordial origin would explain the 
different locations of the gas and the dust. Based on our ALMA 
observations, we propose that HD 21997 is a hybrid system, 
where primordial gas is accompanied by secondary debris dust. 
This fact challenges the current paradigm of disk evolution, 
because the age of HD 2 1 997 exceeds both the model predictions 
for disk clearing and the ages of the oldest T Tauri-like or 
transitional gas disks in the literature (Kastner et al. 2008). 
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